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In SU(3) lattice QCD formalism, we propose a method to extract gauge fields from link-variables analytically. With this method, 
we first study effective mass generation of off-diagonal gluons and infrared Abelian dominance in the maximally Abelian (MA) 
gauge in the SU(3) case. Using SU(3) lattice QCD, we investigate the propagator and the effective mass of the gluon fields in the 
MA gauge with U(l)3XU(l)s Landau gauge fixing. The Monte Carlo simulation is performed on 16 4 at /3 — 5.7, 5.8 and 6.0 at the 
quenched level. The off-diagonal gluons behave as massive vector bosons with the approximate effective mass M Q g =s 1.0-1 .2GeV 
in the region of r — 0.2 - 0.8fm, and the propagation is limited within a short range, while the propagation of diagonal gluons 
remains even in a large range. In this way, infrared Abelian dominance is shown in terms of short-range propagation of off-diagonal 
gluons. Furthermore, we investigate the functional form of the off-diagonal gluon propagator. The functional form is well described 
by the four-dimensional Euclidean Yukawa-type function e~ nkmr jr with m s - 1.5 — 1.6GeV for r = 0.1-0.8 fm. This also indicates 
that the spectral function of off-diagonal gluons has the negative-value region. 
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1. Introduction 

Quantum chromodynamics (QCD) is the fundamental gauge 
theory of the strong interaction based on quarks and gluons. 
There are a variety of nonperturbative phenomena in low energy 
QCD such as color confinement and chiral symmetry breaking. 
These nonperturbative phenomena have been studied both in 
analytical frameworks and in lattice QCD 

CHI- 

On the quark-confinement mechanism, Nambu, 't Hooft 
and Mandelstam suggested the dual-superconductor picture JU]. 
This picture is based on the electromagnetic duality and the 
analogy with the one-dimensional squeezing of the magnetic 
flux in the type-II superconductor. In this picture, there occurs 
color magnetic monopole condensation, and then the color- 
electric flux between the quark and the antiquark is squeezed as 
a one-dimensional tube due to the dual Higgs mechanism. From 
the viewpoint of the dual-superconductor picture in QCD, how- 
ever, there are two assumptions of Abelian dominance |5l |6J] 
and monopole condensation. Here, Abelian dominance means 
that only the diagonal gluon component plays the dominant role 
for the nonperturbative QCD phenomena like confinement. 

The maximally abelian (MA) gauge has mainly been inves- 
tigated from the viewpoint of the dual-superconductor picture 
B7I 417I1 and the various lattice QCD Monte Carlo simulations 
show that the MA gauge fixing seems to support these assump- 
tions uni. 
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According to these studies, the diagonal gluons seem to 
be significant to the infrared QCD physics, which is called 
"infrared Abelian dominance". Infrared Abelian dominance 
means that off-diagonal gluons do not contribute to infrared 
QCD. Therefore, the essence of infrared Abelian dominance is 
the behavior of the off-diagonal gluon propagator. 

The gluon propagators in the MA gauge has been investi- 
gated in SU(2) lattice Monte Carlo simulations (30. To 
investigate the gluon propagators in the MA gauge, the glu- 
ons would be extracted exactly from the link-variables, because 
the link-variable cannot be expanded for a small lattice spac- 
ing due to the large fluctuation of gluons. In SU(2) lattice 
case, the extraction is easy to be done without any approxima- 
tion, because of the SU(2) property. With this extraction, the 
SU(2) lattice simulation suggests that the off-diagonal gluons 
do not propagate in the infrared region due to the effective mass 
M s — 1.2GeV, while the diagonal gluon widely propagates 

m. 

The aim of this paper is to propose a method to extract the 
gluons from the link-variable directly and generally in SU(3) 
lattice QCD, and to investigate the gluon propagators in the MA 
gauge. 



2. Formalism to extract gluon fields from link- variables 

In the Landau gauge, the link-variables can be expanded for 
the small lattice spacing within the scaling region and the glu- 
ons are easily obtained from link-variables approximately 118 - 
2 111 . In Euclidean QCD, the Landau gauge is defined by mini- 
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mizing the global quantity, 

R = J d 4 xA;(x)A a M (x), 



(1) 



by the gauge transformation, where A^(x) e R (a — 1 ■ ■ ■ 8) are 
gluon fields. This means that the gauge fluctuation is maximally 
suppressed in the Landau gauge. On the lattice, the gauge fields 
are expressed by the link-variables U^x) = e mgAf,< - x) eSU(3) 
with the lattice spacing a and the gauge coupling constant g. 
Therefore, in the Landau gauge, \agAn(x)\ <K 1 would be sat- 
isfied within the scaling region and the gluon fields would be 
extracted from the link-variables, 



1 



A„(x) a — [f/ (x) - Ul(x)] - (trace part). 
2iag 



(2) 



In other gauge, however, it is not straightforward to extract the 
gluons from the link-variables. 

In this section, we consider a useful and general method to 
extract the gauge fields analytically and exactly from the link- 
variables whether lagA^x)] <K 1 is satisfied or not lfl9ll . 

To this end, we first define the hermite matrix, 



A . l(„-tf.) 



-{e iagA -e- iagA ) = smagA. 
2i 



(3) 



For simplicity, we have omitted the Lorentz index and space- 
time arguments. 

Arbitrary hermite matrix A can be diagonalized by a unitary 
transformation, 



A,, = QAQ T = 



( M 







A 2 







h ) 



(4) 



where Q eSU(3). We can obtain the eigenvalues A, (i = 1, 2, 3) 
by solving 



det(xl - A) = 0. 



(5) 



This is a cubic equation on x. The eigenvalues A,, i.e, The solu- 
tions of the equation are 



xo,± = zo,± 
where 

= e m 



■ a, 



zo = 



+ e 



-S9/3 



z = e K9±2n)/3 



e 



-i(9±2ff)/3 



(6) 



(7) 



The derivation of the solution and the notation for Eqs. © and 
(Q are given in Appendix A Thus, A t is obtained. 

The unitary matrix O can be also derived as follows. By 
solving Met - A$i, we obtain eigenvectors e; = \xi,yi,Zi) (i = 
1,2,3), |e| = 1. We assume that Zi is nonzero without loss of 
generality and rescale it by 1 \zu 



(8) 





' xdzi ' 
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ydzi 


= A t 


yt/zi 




{ i J 




{ i J 



This is solved easily as 



Xi/Zi 

yt/zi 



-{(An - Ai)(A 22 - Ai) - A 12 A 2 i 



A22 - h -A12 
-A 2 i A u -Ai 



An 

A23 



(9) 



From the normalization condition |et] = 1, we obtain O 1 = 
(el, el, el)- 

When we diagonalize M with the unitary matrix Q, A is also 
diagonalized, 



Aw 



QAO 1 



sin (agQAQ ' ) = 



A x 

A 3 ) 

( sin 6<i 



sin 9 2 







sin 63 t 



(10) 



where -tt/2 < < n/2 (i = 1, 2, 3) is taken. Therefore we can 
derive gluon fields A from link-variables U analytically, 

0i ^ 

3 
( 01 



QAQ. 



1 

ag 



1 

ag 





02 

03 ) 



n. 



(ii) 



This formalism is quite general, because the derivation is cor- 
rect with any gauge and even without any gauge fixing. 

3. SU(3) lattice QCD results for gluon propagators in the 
MA gauge 

Using the SU(3) lattice QCD, we calculate the gluon propa- 
gators in the MA gauge with the U(l)3xU(l)g Landau gauge 
fixing. In the MA gauge, to investigate the gluon propaga- 
tors, we use the gluon fields extracted directly from the link- 
variables. The Monte Carlo simulation is performed on the 16 4 
lattice with [} = 5.7, 5.8 and 6.0 at the quenched level. All mea- 
surements are done every 500 sweeps after a thermalization of 
10,000 sweeps using the pseudo heat-bath algorithm. We pre- 
pare 50 gauge configurations for the calculation at each /3. Error 
is estimated with the jackknife analysis. The MA gauge fixing 
is performed by the maximization of 



Rma = 



(12) 



where H = (T3, T&) is the Cartan generator. In this gauge fix- 
ing, there remains U(1)3XU(1)8 gauge symmetry. After the 
Cartan decomposition for the SU(3) link-variables as t/ M (x) = 
M^u^x) with u^x) = e '(» ! W 3 +» 8 W^) eU(l) 3 xU(l) g and 
Mfi(x) = e ^«*3, 8 0W eSU(3)/U(l) 3 xU(l) 8 , the residual gauge 
fixing is performed by the maximization of 



KucijL = £j X Re ^V*)]- 



(13) 



At j6 = 5.7,5.8 and 6.0, the lattice spacings a are estimated as 
a - 0.186fm, 0.152fm and 0.104fm, respectively, which lead to 
the string tension cr - 0.89GeV/fm in the inter-quark potential 

Here, we study the Euclidean scalar-type propagator of the 
diagonal (Abelian) gluon as 



0=3,8 



and that of the off-diagonal gluon as 



(14) 



(15) 



o#3,8 



These scalar-type propagators are expressed as the function of 
the four-dimensional Euclidean distance r = yOfy - y^) 2 - 

We show in Fig. [T]the lattice QCD result for the diagonal 
gluon propagator G^ el (r) and the off-diagonal gluon propa- 



gator G°"(r) in the MA gauge with the U(1) 3 XU(1) 8 Landau 
gauge fixing. In the MA gauge, G^ el (r) an d G°^(r) mani- 



festly differ. The diagonal-gluon propagator G^j (r) takes a 
large value even at the long distance. In fact, the diagonal glu- 
ons A 3 , A® in the MA gauge propagate over the long distance. 
On the other hand, the off-diagonal gluon propagator G°JJ(r) 
rapidly decreases and is negligible for r Z 0.4 fm in comparison 
with G^ el (r)- Then, the off-diagonal gluons A£ (a + 3, 8) seem 
to propagate only within the short range as r t 0.4 fm. Thus, 
"infrared abelian dominance" is found in the MA gauge. 
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Figure 1: The SU(3) lattice QCD results of the scalar-type gluon propaga- 
tors G^ d (r) and G°Jf(r) as the function of r = J(Xp - y^) 2 in the MA gauge 
with the U(1)3XU(1)8 Landau gauge fixing in the physical unit. The Monte 
Carlo simulation is performed on the 16 4 lattice with fl = 5.7,5.8 and 6.0. 
The diagonal-gluon propagator G^J (r) takes a large value even at the long 
distance. On the other hand, the off-diagonal gluon propagator G°J(r) rapidly 
decreases. 



4. Estimation of the off-diagonal gluon mass in the MA 
gauge 

Next, we investigate the effective gluon mass. We start from 
the Lagrangian of the free massive vector field A M with the mass 



M + in the Proca formalism, 

I = ^Ay-dyA^f+^A^, (16) 

in the Euclidean metric. The scalar-type propagator G / , A ,(r; M) 
can be expressed with the modified Bessel function K\{z) as 



G m {r,M) = (A^A^y)} 
d 4 k 



-! 



= 3 



> 4 k 



,ik-(x-y) 



1 



^_ 

k 2 +M 2 \ + M 2 



f d4k JHx-y) \ , _Lg*( x _ y) 

J {2nf k 2 + M 2 M 2 ° ( y) 

iflT Ki(Mr)+ ^ (x - y) - (17) 



M 2 

In the infrared region with large Mr, Eq. ( TT7b reduces to 

3 VM e- Mr 

G m (r;M) - - — , 

2(27r)i A 

using the asymptotic expansion, 



K 



^ £(|+n) l 

2z 6 ^n\T(l-n)(2z)"' 



i oo 



(18) 



(19) 



for large Re z. 
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Figure 2: The logarithmic plot of r 3 ' 2 G°J(r) and r 3 ' 2 G^ c '(r) as the function 
of the distance r in the MA gauge with the U(1)3XU(1)8 Landau gauge fixing, 
using the SU(3) lattice QCD with 16 4 at/3 = 5.7, 5.8 and 6.0. 



In Fig. |2l we show the logarithmic plot of r 3 ^ 2 G°JJ(r) and 
r'' 2 G^\r) as the function of the distance r in the MA gauge 
with the U(l)3XU(l)g Landau gauge fixing. From the lin- 
ear slope on r 3 / 2 G°JJ(r), the effective off-diagonal gluon mass 
M Q ff is estimated. We summarize in Table 1 the effective off- 
diagonal gluon mass M en - obtained from the slope analysis at 
P = 5.7,5.8 and 6.0. Therefore, the off-diagonal gluons seem 
to have a large mass M a g - 1.0-1 .2GeV. This result approxi- 
mately coincides with SU(2) lattice calculation lfl6Tl . 
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Table 1: Summary table of conditions and results in SU(3) lattice QCD. In the 
MA gauge, the off-diagonal gluons seem to have a large effective mass M g =* 
1.0— 1.2GeV and the functional form in the range of r = 0.1 - 0.8fm is well de- 
scribed with the four-dimensional Euclidean Yukawa function ~ exp(-m ff r)Jr 
with m„(f = 1.5— 1.6GeV. 



lattice size 




o[fm] 


M off [GeV] 


w ff[GeV] 




5.7 


0.186 


1.0 


1.5 


16 4 


5.8 


0.152 


1.0 


1.5 




6.0 


0.104 


1.2 


1.6 



Finally in this section, we discuss the relation between in- 
frared abelian dominance and the off-diagonal gluon mass. Due 
to the large effective mass M g, the off-diagonal gluon propaga- 
tion is restricted within about M~i 0.2 fm in the MA gauge. 
Therefore, at the infrared scale as r » 0.2 fm, the off-diagonal 
gluons (a + 3, 8) cannot mediate the long-range force, and 
only the diagonal gluons A?, A® can mediate the long-range 
interaction in the MA gauge. In fact, in the MA gauge, the 
off-diagonal gluons are expected to be inactive due to the large 
mass M fF in the infrared region in comparison with the diag- 
onal gluons. Then, infrared abelian dominance holds within 
r » M~g. 

5. The analysis of the functional form of the off-diagonal 
gluon propagator in the MA gauge 

In this section, we investigate the functional form of the off- 
diagonal gluon propagator in the MA gauge in SU(3) lattice 
QCD. In the previous section, we compare the gluon propa- 
gator with the massive vector boson propagator and estimate 
the gluon mass. In fact, the gluon propagator would not be 
described by a simple massive propagator Eq. (TTTb in whole 
region of r — 0.1 - 0.8fm. 



There is the similar situation in the Landau gauge 12 111 . The 
functional form of the gluon propagator cannot be described 
by exp(-Mr)/r^ 2 with an effective mass M in whole region of 
r — 0.1 - l.Ofm. The appropriate form is the four-dimensional 
Euclidean Yukawa-type function exp(— mr)/r with a mass pa- 
rameter m. 

In the same way, in the MA gauge, we also compare the 
gluon propagator with the four-dimensional Euclidean Yukawa 
function. In Fig. [3] we show the logarithmic plot of rG°JJ(r) 
and rG^ d (r) as the function of the distance r in the MA gauge 
with the U(1)3XU(1)8 Landau gauge fixing. Note that the 
logarithmic plot of rG°jJ(r) are almost linear in the region of 
r - 0.1 - 0.8fm and then the off-diagonal gluon propagator 
is well expressed by the four-dimensional Euclidean Yukawa 
function in this region, 



(20) 



with a mass parameter m ff and a dimensionless constant A. The 
best-fit mass parameter m g is given in Table 1 at each /3 = 
5.7,5.8 and 6.0. 




Figure 3: The logarithmic plot of rG° 



If) as the function of 



the distance r in the MA gauge with the U(1)3XU(1)8 Landau gauge fixing, 
using the SU(3) lattice QCD with 16 4 at /? = 5.7, 5.8 and 6.0. For rG°£(r), the 
approximate linear correspondence is found. 



We comment on the four-dimensional Euclidean Yukawa- 



type propagator 112 111 . If the functional form of the off-diagonal 
gluon is well discribed by the four-dimensional Yukawa func- 
tion, we analytically calculate the off-diagonal zero-spatial- 
momentum propagator, 



d 3 xGf u (r), 



(21) 



and obtain the spectral function p(io) by the inverse Laplace 
transformation. Also in the MA gauge, the spectral function is 
found to have the negative-value region as in the Landau gauge 

mm in, 



p(oj) = - 



4nAm fi 

( W 2 - m2 ff )3/2 

47rA/ V2m ff 



6(a> - m oB ) 



(oj - m oS y/ 2 



8(0) - m oB ). 



(22) 



6. Summary and Concluding Remarks 

We have first studied the gluon propagators in the MA gauge 
with the U(1)3XU(1)8 Landau gauge fixing using the SU(3) 
lattice QCD. To investigate the gluon propagators in the MA 
gauge, we have considered to derive the gluon fields from the 
SU(3) link-variables. In this method, the gauge fields have been 
extracted by diagonalizing the link-variables and taking the log- 
arithm. Owing to this method, any quantity expressed by the 
gluon fields can be calculated directly from link-variables, even 
if [agA^x^ <s 1 does not satisfy. As one of the general merits 
of this method, we can directly check the correspondence be- 
tween gluon fields and the continuum gauge fixing in arbitrary 
lattice gauge fixing performed with link-variables. In principle, 
with this method, continuum gauge fixing with gluon fields can 
be also performed directly. 
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With this method, we have measured the Euclidean scalar- 
type propagators G^{r) of the diagonal and the off-diagonal 
gluons, and found the infrared Abelian dominance. The Monte 
Carlo simulation is performed on the 16 4 lattice with [5 = 
5.7, 5.8 and 6.0 at the quenched level. We have found that the 
off-diagonal gluons behave as massive vector bosons with the 
effective mass M off 1.0-1.2 GeV for r = 0.2 - 0.8 fm. 
The effective gluon mass has been estimated from the linear 
fitting analysis of the logarithmic plot of r 3 ^ 2 G°^(r). Due to 
the large value, the finite-size effect for the off-diagonal gluon 
mass is expected to be ignored. The large gluon mass shows 
that the off-diagonal gluons cannot mediate the interaction over 
the large distance as r » M~^, and such an infrared inactivity of 
the off-diagonal gluons would lead infrared Abelian dominance 
in the MA gauge. 

On the other hand, from the behavior of the diagonal gluon 
propagator G^ el (r) an d i 3 G^~(f), the diagonal gluons seem 
to behave as light or massless particles. However, for the de- 
tailed argument on Gjr el (r), one should consider the finite size 
effect more carefully, because the diagonal gluons would prop- 
agate over the long distance beyond the lattice size. 

Finally, we have also investigated the functional form of the 
scalar-type propagator in the MA gauge. We show that the 
off-diagonal gluon propagator is well described by the four- 
dimensional Euclidean Yukawa-type form with the mass pa- 
rameter m a ff =i 1.5-1.6 GeV in the region of r = 0.1 - 0.8 
fm. Then, this indicates that the spectral function p(u>) of the 
off-diagonal gluons in the MA gauge has the negative-value re- 
gion as in the Landau gauge. 

On the other hand, the functional form of the diagonal gluon 
propagator seems to be the four-dimensional Euclidean Yukawa 
function with the lighter mass parameter. However, to discuss 
the functional form clearly, the finite size effect is to be checked 
carefully just like the estimation of the diagonal effective gluon 
mass. 

In this study, we investigate the scalar-type off-diagonal 
gluon propagator. To be strict, the off-diagonal gluon propaga- 
tor consists of two scalar functions corresponding to longitudi- 
nal and transverse components. Therefore, we will investigate 
each effective mass and the functional form of these compo- 
nents. 
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Appendix A. The eigenvalues of 3x3 hermite matrix 

In this Appendix we derive the solution for the cubic equa- 
tion, 



where A is an hermite matrix. 





x- An 


-A12 


-A13 


det(xl - A) = 


-A 2 i 


x - A 22 -A 23 




-A 3 i 


-A32 


x- A 33 




x 3 - x 2 (A 


1 1 + A22 


+ A 33 ) 



+ x(A 22 A33 + A33A11 + A11A22 

- |A23| 2 -|A 3 i| 2 -|A 12 | 2 ) 

- detA = 

=> x 3 - 3ax 2 -fix - y = 0, 



(A.2) 



where 3a = TrA e R, /3 = -(A22A33 + A33A11 + A11A22 - 
|A 23 | 2 - IA31I 2 - |A 12 | 2 ) e R, y = detA e R. Here, when we 
define x — y + a, this equation is rewritten, 

y 3 - (3a 2 + p)y - (2a 3 + a/3 + y) = 0, 



where 



3a 2 +y = -{(A n 
6 



(A.3) 



■ A 22 ) 2 + (A 22 - A33) 2 + (A33 - A„) 2 
+6(|A 2 3| 2 + |A3i| 2 + |A 12 | 2 )}. 



For simplicity, we define p = yja 2 + /3/3 e R, q = 2a 3 + a/3 + 
ye R, 

y 3 -3p 2 y-q = 0. (A.4) 

If p + 0, we can rescale it by 1 / p. When we define y — pz, this 
equation is rewritten, 

z 3 - 3z - 2r = 0, 



where 2r = qp 3 e R . When z — u + v, 

u 3 + v 3 + 3uv(u + v) - 3(u + v) - 2r = 
<=> u 3 + v 3 - 2r + 3(u + v)(uv - 1) = 0. 



(A.5) 

(A.6) 
(A.7) 



Without loss of generality, we can determine uv — 1 and obtain 



w t v — 2r 



u\ 3 



I. 



These solution are obtained by solving 

f 2 - 2rt + 1 = 0. 
Therefore this solution is 
t 



= r ± i VI 

- „±<0 



(A.8) 
(A.9) 

(A. 10) 
(A.ll) 



where r 2 < 1 could be proven. By using this, z can be obtained. 
Note that the solutions of z are real, because the eigenvalues of 
A are real. 



-a/3 



det(xl - A) = 0, 



(A.l) 



zo = er- +e ~- (A. 12) 

z± s e *fe*)/3 +e -*fea0/3. (A.13) 

Therefore we can obtain the solutions of Eq. (IA.U . 

xo,± = pzo,± + a 

= zo,± ^a 2 +/3/3+a. (A. 14) 
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